Finite-Size effects in "Single Chain Magnets": an experimental and theoretical study 
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The problem of finite size effects in s=l/2 Ising systems showing slow dynamics of the mag- 
netization is investigated introducing diamagnetic impurities in a Co^^-radical chain. The static 
magnetic properties have been measured and analyzed considering the peculiarities induced by the 
ferrimagnetic character of the compound. The dynamic susceptibility shows that an Arrhenius law 
is observed with the same energy barrier for the pure and the doped compounds while the prefactor 
decreases, as theoretically predicted. Multiple spins reversal has also been investigated. 
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One of the open problems in low dimensional mag- 
netism is the influence of finite-size effects 0, 12 and in- 
vestigation of their role on dynamic properties, which are 
of great relevance for nanostructures, is particularly in- 
teresting. For the Ising spin chain Glauber proposed a 
stochastic dynamics 3] widely applied to a variety of phe- 
nomena. The experimental investigation of this dynamics 
in magnetic ID compounds was prevented by the strict 
requirements the system must fulfill. Recently molecu- 
lar ID materials were found to show magnetic hysteresis 
in absence of long range magnetic order 0, 0, They 
have thus been named "Single Chain Magnets" (SCMs) 
by analogy to slow relaxing OD systems widely known 
as "Single Molecule Magnets" (SMMs)[3. The first of 
these compounds was Co(hfac)2NitPhOMe (CoPhOMe 
in the following) "t] , made of alternating and interact- 
ing ions (anisotropic effective s=l/2 spins) and 
nitronyl-nitroxide radicals (isotropic s=l/2 spins). The 
relaxation time follows an Arrhenius law as predicted by 
Glauber 3] for a ID Ising system, with r = r^e^/^^'^ and 
A/Zcs =152K, for ten decades of time in the temperature 
range T =4-15K [|. 

The observation of Glauber dynamics in the paramag- 
netic phase of a real system represents a paradox. In 
fact a strong Ising type intra-chain interaction is re- 
quired because A is proportional to the coupling con- 
stant J. On the other hand the correlation length of the 
Ising model diverges exponentially at low temperature, 
^ ~ exp{j^), and favors the 3D ordering by enhancing 
the weak inter-chain interactions |^ . A possible explana- 
tion could be the presence of naturally occurring defects, 
which can geometrically limit the correlation length and 
dramatically reduce the 3D ordering temperature 0, llfl| | . 
The presence of defects is known to affect also the dy- 
namics of the magnetization, but direct investigations of 
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their influence on the slow relaxing magnetization are 
still lacking. The persistence of slow dynamics even for 
small segments of chains is an important result, mak- 
ing SCMs much more appealing than SMMs in many 
respects. Finite-size effects are commonly investigated 
through the insertion of breaks by doping with diamag- 
netic ions, but important aspects, like the reproducibility 
or homogeneity of the concentration of the dopant, are 
often overlooked. We present here a thorough experimen- 
tal and theoretical investigation of finite-size effects in the 
slow relaxing Ising chain CoPhOMe, thus evidencing the 
key role played by the defective sites in the nucleation of 
the excitation, similarly to what is observed in quasi-lD 
oxides 0. A short-cut of the Glauber barrier through 
coherent reversal of multiple spins is observable in the ac 
susceptibility thanks to the introduction of defects. 

Variable concentrations of diamagnetic impurities have 
been inserted in the chains and the samples have been in- 
vestigated with unprecedented accuracy. Doped crystals 
were prepared [ij] adding to the solutions different ra- 
tios of Zn(hfac)2 • 2H2O vs Co(hfac)2 -21120 in the range 
0-0.30. As a strongly inhomogeneous distribution of the 
dopant inside single crystals could be of major relevance 
for the magnetic properties, we used the Particle In- 
duced X-ray Emission technique (PIXE) with an external 
micro-beam (12j| to analyze the Zn doping profile. The 
low minimum detection limits of PIXE allowed us to find 
that no metals other than Co and Zn were present in con- 
centrations >12 ppm and thus no appreciable contribu- 
tion could rise from paramagnetic impurities. The mean 
Zn/Co ratios of the investigated samples, c, were 0.3%, 
1.9%, and 4.7%, giving an average length L between two 
Zn^+ ions variable from ca. 300 to ca. 20 spins. The Zn 
concentration was always found to be lower than the ratio 
of the starting solution and reproducible in crystals from 
the same batch (within an error of 5% on c) but not repro- 
ducible in different syntheses. Using the micro-beam we 
could then monitor the spatial distribution of the dopant 
inside individual crystals, a procedure never performed 
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FIG. 1; (Color online) TjajCo ratio in a single crystal of 
CoPhOMe. The distribution of the dopant is reported 
for both longitudinal and transversal directions, scanned as 
shown in the drawing. Lines are guides to the eye. 

before in this kind of studies. The resuhs are shown in 
Fig.l for longitudinal and transversal scans. The former 
shows uniformly distributed impurities while the latter 
shows that the Zn concentration slightly increases going 
from the center to the edges. This observation is con- 
sistent with a progressive enriching of the solution in Zn 
content. 

The magnetic susceptibility of a chain in a moderate 
field is strongly sensitive to geometrical limitations of the 
correlation length. In Fig. 2 (top) we report the tem- 
perature dependence of the real component of the ac 
magnetic susceptibility measured at 2.7kHz in a static 
magnetic field H = 2kOe along the easy axis for an un- 
doped crystal and for two doped samples. In contrast 
with magnetic data and theoretical calculations reported 
so far for cha ins, which show the presence of only one 
peak El El mi, it is evident the presence of a large 
anomalous structure with two peaks for T ^ 33K and 
for T ~ 14K (below the blocking temperature Tb ~ 12K 
we observe frequency dependent dynamical effects). The 
higher temperature peak roughly corresponds to that ex- 
pected for an infinite ID Ising system in presence of a field 
(see Fig. 2 bottom) and shifts to higher temperature on 
increasing the field, as expected. This peak is more pro- 
nounced in the pure sample and decreases and shifts to 
lower temperatures with higher dopings, while the lower 
temperature peak increases. At high doping (0.047) only 
this latter peak survives. This unprecedented behavior 
can be explained solving the randomly diluted ferrimag- 
netic Ising chain in field. The spin Hamiltonian we have 
used is: 

= - ^[J(cr2i_icr2i+0'2iCr2i+l)H ^ — (5CoO'2i+g_RCr2i-l)] 

i=l 

(1) 

where gco and qr are the Lande factors of Co and radical 
spins, /is is the Bohr magneton, H is the external mag- 
netic field and <Ti = ±1, except for (T^+i = 0. Within a 
transfer matrix formalism, which is reported in literature 
[T^ Q, ^3 , the free energy of a segment of length L can 




10 20 30 40 50 

T(K) 



FIG. 2: (Color online) (top) Temperature dependence of the 
longitudinal magnetic ac susceptibility (2.7kHz) measured for 
different dopings in a static field of 2kOe with an ac field of 
10 Oe. Data for other frequencies in the 10 Hz-lOkHz range 
are not reported because x' is frequency dependent only in 
the shaded area, (bottom) Calculated susceptibility (param- 
eters in the text) for different concentrations of impurities. A 
parabolic distribution is used, at low concentrations, to ac- 
count for the observed transversal trend. 



be obtained. The thermodynamic observables are calcu- 
lated summing all the Us weighted by the probability of 
occurrence of that length: = c^(l — c)^ [T^ . 

In Fig. 2 (bottom) we report the calculated magnetic 
susceptibility x{T) for J/fcs=-90K, iJ=2kOe, gco^^ and 
gR=2 and for different doping values. For the infinite 
chain we obtain only one peak at T ~ 34K, which shows 
the observed trend with field, while the introduction of 
non-magnetic impurities gives rise to the observed peak 
shift and to a shoulder at about 15K. This feature, which 
eventually develops into a single peak at high doping, 
is mainly due to the finite-size contribution to the free 
energy [iJl, as it will be shown in a forthcoming work. 
Although complete quantitative agreement is not possi- 
ble due to the helicoidal structure of the CoPhOMe chain 
1^ ^3 the overall behavior and the anomalous feature 
are well reproduced by the simple model expressed by 
(1). This is especially true if one introduces, at low con- 
centrations, a parabolic distribution of the dopant, as in 
Fig. 1 , while at higher concentrations the behavior can be 
reproduced using either a distribution or homogeneous 
doping. 

The presence of the low temperature peak even in the 
pure compound suggests that other type of defects are in- 
deed present beyond the introduced Zn impurities. Both 
crystalline defects and chemical modifications (due to a 
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FIG. 3: (Color online) Relaxation time of the magnetization 
extracted from the longitudinal ac susceptibility measured in 
zero static field for different concentrations of Zn. 
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FIG. 4: (Color online) Calculated temperature dependence of 
the relaxation time extracted from x('^i T) with ao = 3.6-10^'' 
sT^ and q — 0.73. To compare with data in Fig. 3 assume 
J/ks^m K. 

certain lability of the organic radicals) are present in the 
undoped compound and are expected to be more con- 
centrated close to the crystal surface, as for the dopant. 
These data suggest a concentration of defects in the pure 
compound of the order of a few per thousand, a value 
that agrees well with the calculated curves of Fig. 2 and 
the following discussion on the dynamics. 

The dynamics of the magnetization of the pure and 
doped samples has been investigated measuring the tem- 
perature dependence of the real (%') and imaginary (x") 
components of the ac susceptibility for seven different 
frequencies in the range lOHz-lOkHz. x" displays a max- 
imum whose position depends on the frequency w of the 
applied ac field, and the relaxation time can thus be ex- 
tracted 0: T~^{Tynax) — w. For the highest concen- 
tration the curves become more distorted, with a small 
shoulder at low temperature. The results are plotted in 
Fig. 3 in the ln(T) vs 1/T scale. A linear behavior is ob- 
served for all compounds and the slope {i.e. the barrier 



A of the Arrhenius law) remains substantially unchanged 
(160 ± 8K). The fitted lines shift down on increasing 
the doping with the pre-exponential factor tq going from 
3.5 • 10~^^s for the pure compound to 1.0 • lO^^^s when 
c=0.047. 

In order to reproduce theoretically the data of Fig. 3 
we assumed the simplified model of a dilute ferromagnetic 
chain (in a ferrimagnet only an equivalent ferromagnetic 
branch with J = \ J\ and g = gco — 9r is effective for the 
slow relaxation!^). For an open chain of length L the 
time evolution of the expectation values of the spins can 
be described by the matrix equation: 

|S = -MS (2) 

where M is a tridiagonal matrix. The time scale of the re- 
laxation process (T) is related to the smaller eigenvalue 
of M. In the range of ^ and L ^ 1 asymptotic expan- 
sions for this eigenvalue are given in literature, leading 
to tl{T) - if C < T, and tl{T) - if ^ > l[i1|U 
However, as the experimental t is extracted from the ac 
susceptibility of an ensemble of segments with different 
length, we need to evaluate tl for any L and T . The 
eigenvalue spectrum of M can be expressed as 

\L{e) = qaa[l - -icos{e)\ (3) 

where q is the probability of reversal of the isolated spin 
per unit time a^^ , 7 = tanh(j^^), and by 9 we denote 
the L possible roots of the trigonometric equation 2.6 of 
ref. (13 • The relaxation time of the segment of length 
L is related to the smallest nonzero root via tl{T) = 
X~[^{9o), while the root 6 — is rejected since it provides 
the relaxation time of the infinite chain. We introduced 
tl{T) in the average susceptibility: 

where xl is the susceptibility of a ferromagnetic Ising 
segment of length L [lj|. We extracted the relaxation 
time of the randomly doped sample from the maximum 
of the imaginary component of the computed ac suscep- 
tibility. In Fig. 4 we show the Arrhenius plots thus ob- 
tained for different concentrations of dopant. A crossover 
between the two regions with L ^ and L <C ^ is ob- 
served and the slope is halved in the latter case, where 
the concentration of impurities c only affects the offset 
of the curves. This is consistent with a relaxation pro- 
cess where a domain wall is nucleated at an end point of 
a segment L, with a halved energy cost 2 J, and has the 
probability to reach the other end point jl^j leading to 
the linear dependence of the relaxation time on the size 
of the system. 

The experimental data of Fig. 3 scale with c as the cal- 
culated curves when L <C ^, thus evidencing that finite 
size effects are also present in the nominally pure com- 
pound. Although a direct analysis of its defects is not 
possible, by comparing the values of tq we estimate an 
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FIG. 5: (Color online) (a) Calculated Cole-Cole plot for dif- 
ferent concentrations of randomly distributed diamagnetic 
impurities at T = 8K assuming J/fcs = 80K, q = 0.73, 
ao = 3.6 • 10" s"\ The half-filled circles stay for c=0.047 
without multiple spin reversal, (b) Experimental data at 
T = 8K for the four investigated samples. 

average length of the chains in the undoped sample of the 
order of a thousand spins (as deduced from the analysis 
of the susceptibility in a moderate field). The geometri- 
cal reduction of ^ can be responsible for the suppression 
of long range magnetic order in CoPhOMe 0, • The 
measured activation energy A thus corresponds to 2 J in- 
stead of 4J, leading to J /ks = 80 K, in good agreement 
with the simulation of the susceptibility data of Fig. 2. 

Up to now we have neglected a relaxation mechanism 
which involves the simultaneous reversal of all the L spins 
of a segment, whose probability scales as and therefore 
is relevant for short segments. This term can be added 
to Xl{0{)) to obtain the probability to relax in the unit 
time Qfn ^- Clear evidences of such a relaxation mecha- 



nism appear in the Cole-Cole plot (Fig. 5a), where x" vs. 
x' is reported. In the absence of a distribution of relax- 
ation times the curve is expected to be a semicircle with 
the center on the x axis. A distribution in t, induced by 
the random doping, has the only effect to slightly push 
the center down, even for the highest doping. On the 
contrary, the coherent reversal of all the spins gives rise 
to a more complex behavior. The calculated curves of 
Fig. 5a are in fact well below the semicircle and show a 
flat region for higher frequencies where x' ^ 0. This 
effect is particularly evident for high doping. The exper- 
imental results are shown in Fig. 5b. Although the range 
of available frequencies is limited, the agreement between 
observed and calculated curves suggests that coherent re- 
versal of all the spins is experimentally observed for short 
segments. An accurate investigation on the temperature 
dependence of q could discriminate between a thermally 
activated and a tunnelling process. 

In conclusion we have shown that finite size effects are 
important in magnetic bistable ID nanostructures halv- 
ing the energy barrier predicted by the Glauber model. 
We observe a linear dependence of the relaxation time 
on the size of the system, which is completely different 
from the exponential dependence on the square of the 
spin value observed in SMMs, or the exponential depen- 
dence on the volume in single domain particles. This 
opens the possibility to employ segments of chains, even 
without a precise control on the length, as nanometric 
magnetic memory units characterized by a higher T;, and 
a larger magnetic moment than those obtained in SMMs. 
Moreover a non thermally activated relaxation could be 
observed for short segments. To better investigate this 
phenomenon higher doping is required but it cannot be 
obtained with the present chemical approach. Synthesis 
in confined media, as well as the organization of segments 
on surfaces, as already done for metallic Co nanowires on 
Pt surfaces , are under investigation. 

Financial support from Italian MIUR (Fondi FISR and 
FIRB), German DFG, and Brazihan CAPES, CNPq and 
Inst, de Nanociencias is gratefully acknowledged. 



[1] G. Xu et ai, Science 289, 4192 (2000). 
[2] S. Eggert et ai, Phys. Rev. Lett. 89, 047202 (2002). 
[3] R. J. Glauber, J. Math. Phys. 4, 294 (1963). 
[4] A. Caneschi et ai, Angew. Chem., Int. Ed. Eng. 40, 1760 
(2001). 

[5] R. Clerac et al, J. Am. Chem. Soc. 124, 12837 (2002). 
[6] R. Lescouezec et al., Angew. Chem., Int. Ed. Eng. 42, 
1483 (2003). 

[7] D. Gatteschi, R. Sessoli, Angew. Chem., Int. Ed. Eng. 

42, 268 (2003). 
[8] A. Caneschi et al, Europhys. Lett 58, 771 (2002). 
[9] D. Hone et al, Phys. Rev. B, 12, 5141 (1975). 
[10] C. Dupas et al, Phys. Rev. B, 25, 3261 (1982). 



[11] A. Caneschi et al, Dalton Trans., 21, 3907 (2000). 
[12] M. Massi et al, Nucl Instrum. Methods B, 190, 276 
(2002). 

[13] M. Wortis, Phys. Rev. B, 10, 4665 (1974). 
[14] F. Matsubara et al. Can. J. Phys., 51, 1053 (1973). 
[15] M. G. Pini, A. Rettori Phys. Lett. A, 127, 70 (1988). 
[16] J. Kamphorst Leal da Silva et al, Phys. Rev. E, 52, 4527 
(1995). 

[17] J.H. Luscombe et al, Phys. Rev. E, 53, 5852 (1996). 
[18] R. Cordery et al, Phys. Rev. B, 24, 5402 (1981). 
[19] P. Gambardella et al. Nature 416, 301 (2002). 



